No indications of metal-insulator transition for systems of interacting electrons 

in two dimensions 
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The influence of Coulomb interaction on transport properties of spinless electrons in small disordered two 
dimensional systems is studied within a tight binding model. Spatial correlations, inverse participation ratio, and 
multifractal spectrum of the zero temperature local tunneling amplitude as well as the DC Kubo conductance are 
traced as function of the interaction strength U. When U is increased, all of the above quantities are shifted rather 
smoothly towards localized behavior, indicating the absence of an interaction driven insulator-metal transition. 



PACS numbers: 71.30.+h, 64.60.Ak, 73.20.Fz 

Recently, much interest is focused on the influence of 
electron-electron interaction (eei) on the localization proper- 
ties of two dimensional electrons in disordered systems. It is 

, motivated by the experimental observations of a crossover in 
the behavior of the conductance of low density two dimen- 
sional electrons from an insulating like temperature depen- 
dence at low densities to a metallic one at higher densities [[!]]. 
In some cases a transition back to an insulating dependence 
at even higher densities was observed [^]. This so dubbed 
2DMIT (two dimensional "metal-insulator" transition) is at 
odd with the scaling theory of localization [|J which, for non- 
interacting electrons, asserts that all states in 2D are localized 
by any amount of disorder. It prompted an intensive theoreti- 
cal effort including analytical |Q| and numerical j^-|l0|] work 
which tried to explain the 2DMIT as a result of eei, not taken 
into account by the scaling theory. Another viewpoint, sup- 

■ ported by some recent experiments [pDj, maintains that there 
is no 2DMIT (that is, the systems are insulating at zero tem- 
perature) and that the origin of the observed temperature de- 
pendence is not a result of a metallic zero temperature phase 

& 

In the absence of eei, there are well established relations be- 
tween the conductance and the properties of individual single 
electron states. Such relations do not exist for the many elec- 
tron states. For example, the sensitivity of a single electron 
state at the Fermi energy to a change in the boundary condi- 
tion is a measure of the conductance of the system (the Thou- 
less criteria Jl3|]). Formulation of a similar conjecture for 
many-electron states is much more intricate Jl4| , p3| ]. Another 
example is the behavior of the level statistics at the Ander- 
son localization transition. This transition is accompanied by 
a transition in the single electron energy level statistics from 
Wigner to Poisson [Q, while the many particle energy level 
statistics for non-interacting electrons will remain Poissonian 
both in the metallic and in the insulator regime [ |l7| | . 

In previous studies of 2D spinless electrons with eei [j5j-||] 
numerous equilibrium properties such as many particle energy 
level statistics, two-electron localization length high above the 



Fermi energy, persistent current flow patterns, and charge den- 
sity response to an external perturbation were studied. As 
function of the interaction strength U, the footprints of some 
of these quantities undergo qualitative changes. It might then 
be tempting to interpret it as an evidence for a 2DMIT. How- 
ever, the properties listed above are not directly related to zero 
temperature transport properties of the electron system, which 
are naturally measured in experiments. It is therefore manda- 
tory to carefully examine some quantities which are directly 
related to transport properties of the system. Since the prob- 
lem under study involves disorder and eei the region of in- 
terest in parameter space is that for which both of them are 
significant. This excludes the possibility of using perturbation 
theory, and leaves us with the necessity of employing exact 
numerical diagonalization. 

This task, which is feasible for relatively small systems, is 
carried out in the present Letter. The statistical properties of 
tunneling amplitudes and the conductance distribution of the 
system are examined. It is shown that eei systematically at- 
tenuates the transport through the system and enhances its in- 
sulating features. Thus, although eei may significantly alter 
some of the properties of a 2D electron system as demon- 
strated in previous studies t^)-^, it is conjectured here that 
there is no numerical evidence that it can drive a transition in 
the transport properties of spinless electrons. 

We consider systems composed of 4, 6, and 8 interacting 
electrons residing on 6 x 6, 5 x 5, and 5x4 lattices (cor- 
respondingly) having a torus geometry. It is of course not 
possible here to directly mimic the experimental procedure in 
which the conjectured 2DMIT is driven through variation of 
electron density. Instead, the physical content of this density 
variation can be captured by controlling the ratio of the Fermi 
energy to the interaction energy. In the present model, it is 
achieved simply by changing the interaction strength U while 
keeping other parameters intact [^J. 

In the Tight-Binding approximation, the Hamiltonian of the 
system is given by: 
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where r = (k,j) denotes a lattice site, 
creation operator, ekj is the site energy, chosen randomly be- 
tween — W/2 and W/2 with uniform probability, V is a con- 
stant hopping matrix element and s is the lattice constant. Us- 
ing the Lanczos method we obtain the many-particle eigen- 
values and eigenfunctions \a N ), where N is the number 
of electrons. The zero temperature local tunneling amplitude 
(O^latlO^ -1 ) between the ground state of N and N- 1 elec- 
trons can be employed here in order to characterize the trans- 
port properties of the many-particle interacting system. It has 
the advantage that only the ground state energy and eigenvec- 
tor for N and N — 1 electrons need to be calculated. The 
use of the tunneling amplitude in this context can be moti- 
vated and substantiated by the following considerations: In 
the independent particle approximation, the tunneling density 
of state (TDOS) is given by 
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where if> n is the n-th single particle eigenvector and e„ is the 
n-th single electron eigenvalue. For the many-body interact- 
ing system, the TDOS is defined as [|l8|| : 

"(*) = E K^lE^O^ 1 )! 2 S(e - (e N a - e^ 1 )). (3) 
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The conductance a(e) is related to the transmission t(r, f ' , e) 
of an electron with energy e between two points r, r 1 on 
the interface of the system with external leads through the 
Landauer formula ^ er(e) = (e 2 /fc)^. f , |t(r, r e)| 2 , 
where the sum is over all points on the interface. For a non- 
interacting system the transmission is expressible in terms of 
single particle wave functions, 



(4) 



As discussed in Ref. [pOj], an appropriate Landauer formula 
connecting transmission and conductance for an interacting 
system coupled to external leads has a similar structure, in 
which the transmission in the interacting region is given by 
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In the absence of interaction, Eq. (^|) is reduced to Eq. (Q), 
and Eq. (^]) becomes identical to Eq. (Q). Indeed, the tun- 
neling amplitude (a^latjO^ -1 ) couples only many-particle 
states different from each other by the addition of a single par- 
ticle to an unoccupied single particle state. Thus, for a many- 
particle state \a N ) which corresponds to a Slater determinants 



in which all N — 1 single particle state are occupied as well as 
the n-th single particle state above the Fermi energy, one ob- 
tains ip n (r) = (a^latlO^" 1 ) ande^-e^" 1 = e„, while for 
other many particle states the matrix element vanishes. This 
behavior suggests that the tunneling amplitude is the appro- 
priate quantity to replace the single electron wave function in 
studying transport properties of interacting systems. A similar 
procedure is employed in Ref. [|T]| in order to generalize the 
concept of inverse participation ratio for interacting systems. 

Note, however, that once interactions are present a many 
particle state is a superposition of many different Slater deter- 
minants, and, generically, 5Z f -|(0 JV |a'L|0'' v ' -1 )| 2 1> where 
for states corresponding to quasi-particles the matrix element 
are dominant [g2[|. Thus, for interacting systems the tun- 
neling amplitude (O^latjO^ -1 ) is not normalized. The rea- 
son is, that, due to interaction, the basic objects are quasi- 
particles (rather than particles) with a finite life time. In 
this context, interaction affects transport properties of a given 
system in two different ways. The first (which is not re- 
lated at all to quantum localization) is implied by variation 
of the density of states at a given energy, while the sec- 
ond one (which is the essence of quantum localization) is 
manifested through a change in the correlation properties of 
the tunneling amplitude between different points. Therefore, 
in order to study the influence of eei on quantum localiza- 
tion it is useful to define an effective tunneling amplitude 
m = (O^latlO^-^/CE^IatlO^" 1 } 2 ) 1 / 2 , which will 
be used to study transport properties of an interacting system. 
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FIG. 1.: The decay of the correlation function C(r) for <j>(r) for 
two different disorder strengths W = 5 (in units of V) in the upper 
part and W — 15 in the lower part and 4 interaction strengths (see 
legend). The correlation function decays faster for stronger disorder 
and for stronger interaction. 

The effective tunneling amplitude </>(r) of the system intro- 
duced above is calculated for several values of disorder, and 
then analyzed by employing well known procedures for elu- 
cidating properties of single particle wave functions such as 
spatial correlation function, inverse participation ratio and a 
multifractal analysis. Results for the correlation function 
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for the 6 x 6, 4 electrons system averaged over 100 realizations 
are shown in Fig. [I] (similar results are obtained for the other 
system sizes). Evidently, C(r) decays as function of r indicat- 
ing the loss of amplitude correlations. The decay is faster for 
stronger disorder as is the case for a system of non-interacting 
electrons. Moreover, the decay of C(r) evolves smoothly with 
increasing interaction strength U, indicating stronger localiza- 
tion at higher values of U, showing no sign of a 2DMIT or 
an intermediate metallic phase. These results are significant, 
since an identical system as considered here undergoes a tran- 
sition in the character of its persistent current at intermediate 
values of U [^3jQ] . In particular, the flow pattern of the local 
persistent current is shown to be ordered and essentially one 
dimensional. Our new result indicate that this kind of behav- 
ior should not be interpreted as a signature of a delocalized 
(Coulomb) phase. 

A similar feature is exposed in the inverse participation 
ratio P^ 1 with P = |<^(tD| 4 - Figure ||(a) shows the 
monotonous decrease of P^ 1 with increasing U and disorder 
W. Again, there is no sign of an intermediate metallic phase. 
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FIG. 2.: (a) Participation ratio P _1 as function of interaction 
strength U for 6 x 6, 4 electrons; 5 x 5, 6 electrons; 5 x 4, 8 
electrons; and different disorder values (all averaged over 100 re- 
alizations of disorder). Parts (b) and (c) show the results D(— 2) 
and D(+2) from a multifractal analysis. The results indicate a 
monotonous crossover to stronger localized states with increasing 
interaction or disorder. 

In order to further elucidate the characteristics of the ef- 
fective tunneling amplitude <j>(f) we consider its multifractal 
pattern (for a review, see [|^]). The multifractal exponents 
r(q) are obtained from the scaling behavior of the partition 
function 



(7) 



where the sum runs over all boxes i of linear size £ and pi (£) 
are the probability densities in each box. Figures |^(b) and (c) 
show D(q) = r(q)/(q — 1) for q = —2 and +2 as function of 
interaction strength U. Again, we do not find any indication 
of a 2DMIT. 

Another way to present the results of the multifractal anal- 
ysis is the singularity spectrum f(a), which is the Legendre 



transform of r(q). For each singularity exponent a, f(a) can 
be interpreted as the fractal dimension of the subset of boxes 
i, which are characterized by Pi(£) ~ £ a . Figure [}] shows the 
f(a) spectra for the 6 x 6, 4 electrons system at two disorder 
strengths. Since the singularity spectra are not independent of 
interaction strength (or system size) the multifractality is not 
an indication of a critical point. The spectra become wider 
with increasing interaction or disorder indicating stronger lo- 
calized states, rather than generic multifractal states. Yet, the 
bottom line is that we do not find any indication of a metal- 
lic phase but a mere monotonous crossover into the regime of 
stronger localization. 




2 3 4 

exponent a 



FIG. 3.: Singularity spectra for effective multi particle transmis- 
sion functions (pp. For each singularity exponent a, f(a) can be 
interpreted as the fractal dimension of the subset of boxes i, which 
are characterized by pi (I) ~ t a . The spectra become wider with 
increasing interaction or disorder indicating a crossover to stronger 
localized states. 

In order to relate our results for the tunneling amplitude 
to the standard Kubo formulation of the conductance, we de- 
termine the distribution of the conductances a for different 
disorder configurations and, further, study the dependence of 
the average conductance on the interaction strength U. The 
conductance at zero temperature is calculated using the many- 
particle Kubo formula [|l4j] 
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where 



£ a ~ £ o > Jx is the current operator and M 
is the number of sites. The inelastic broadening 7 is chosen 
to be of the same order as the rescaled "single electron" level 
separation defined as B/N{M — N), where B is the width of 
the many particle energy band [^,^]. The calculation of the 
conductance is cumbersome since it requires the computation 
of many-body eigenfunctions for all the low lying excitations 
(20 in our calculation), which is a difficult numerical task even 
within the Lanczos algorithm. 

We first present the behavior of the conductance as function 
of U. For the 6 x 6, 4 electrons system, the distributions of the 
dimensionless conductance g = (h/e 2 )a are shown in Fig. [| 



for two different values of disorder (W — 5, W = 15). For 
the weaker disorder, which for this small system corresponds 
in the non-interacting case to a marginal metal |{7|], the dis- 
tribution approaches a log-normal distribution as interaction 
strength increases. For the stronger disorder, which even for 
the non-interacting case is in the localized regime, the log- 
normal distribution is preserved for any strength of interac- 
tion. Since the conductance distribution of a localized system 
is expected to follow a log-normal distribution [^6|], while a 
metallic system should follow a normal distribution this be- 
havior indicates that repulsive electron-electron interactions 
suppress the metallic characteristics of the Kubo conductance. 
This is consistent with results reported sometime ago pertain- 
ing to the influence of eei on the conductance [14 2^]. In 
particular it has been shown that in the diffusive regime the 
average conductance is suppressed by interactions, while only 
deep in the localized regime, where the average current is very 
small, some enhancement is possible \ 21\. 
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FIG. 4.: The log of the dimensionless conductance distribution for 
different values of disorder and interaction strength for the 6 x 6, 4 
electrons system. Inset: the average of the log of the dimensionless 
conductance as function of the interaction for the 6 x 6, 4 elec- 
trons system (averaged over 280 realizations of disorder); 5 x 5, 6 
electrons (60 realizations) and 5 x 4, 8 electrons (100 realizations). 

The same trend can be seen when inspecting the average 
logarithmic conductance, displayed in the inset of Fig. 
Thus, neither the average conductance nor its distribution 
show any traces of a 2DMIT. 

In conclusion, we have investigated the transport properties 
of spinless electrons in disordered clusters. Both the tunneling 
amplitude and the conductance do not show any evidence of a 
2DMIT. It has been shown that the tunneling amplitude is the 
natural analog of the single electron wave function appropriate 
for interacting many-particle systems. This is quite attractive, 
since numerous methods developed to connect the properties 



of single electron wave functions to transport properties of the 
system may now be efficiently applied to the tunneling ampli- 
tude function. 
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